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INJECTIVITY AND ALMOST GLOBAL STABILITY OF 
HURWITZ VECTOR FIELDS 

ALVARO CASTANEDA AND VICTOR GUfNEZ 


Abstract. We present, in dimension n > 3, a survey of examples to: the 
Jacobian conjecture, the weak Markus-Yamabe conjecture. Furthermore, we 
show and construct new examples of vector fields where the origin is almost 
globally asymptotically stable by using the novel concept of density functions 
introduced by Rantzer. 


1. Introduction 

One of the central problems on dynamical systems is to determine conditions 
under which certain points or sets are attractors for some dynamics, that is, the 
orbits through points in a neighborhood of the attractor converge to them. In 
the case of continuous-time, that is, flows associated to vector fields, an analytic 
condition ensuring that an equilibrium point x* is a local attractor is given by 
the negativeness of the real part of the eigenvalues of the Jacobian matrix at x*. 
Motivated by this simple observation, in [15], L. Markus and H. Yamabe establish 
their well known global stability conjecture 

Markus-Yamabe Conjecture (MYC): Let F : K." —R" be a vector field 
with F{0) = 0. If for any x G K” all the eigenvalues of JF{x), the Jacobian matrix 
of E at X, have negative real part, then the origin is a global attractor of the system 
X = F{x). 

Let us recall that the vector fields satisfying the hypothesis of MYC are called 
Hurwitz vector fields. It is known that the MYC is true when n < 2 and false 
when n > 3 (see [8] for a counterexample). The proofs in the planar context, both 
the polynomial case (G. Meisters and C. Olech in [16]) as the C^— case (R. Fefiler 
in [m, A.A. Glutsyuk in [13] and C. Gutierrez in [14]) are based on a remarkable 
result of C. Olech [T^, where the author showed that MYC (in dimension two) 
is equivalent to the injectivity of the map F. In R", the problem of knowing if a 
Hurwitz vector field is injective is known as the Weak Markus-Yamabe conjecture. 

Weak Markus-Yamabe Conjecture (WMYC): If F : R" R’" is a 

Hurwitz map, then F is injective. 

The WMYC is true when n < 2 and, to the best of our knowledge, it has been 
proved in dimension n > 3 for Lipschitz Hurwitz maps by A. Fernandes, C. 
Gutierrez and R. Rabanal in [TO] Gorollary 4]. 
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For n = 2, the strong injectivity theorem of Gutierrez [14] is the following: “A 
map / : ^ is injective if [0,oo) fl Spec(J/(a;)) = 0, for all x S R^”. 

However, this result fails in high dimensions as shown B. Smyth and F. Xavier 
EOl Theorem 4]: “There exist integers n > 2 and non-injective polynomial maps 
/ : R" ^ R" with [0, oo) n Spec(J/(a;)) = 0, for all a; G R””. 

In addition, a new tool (density functions) introduced by A. Rantzer in [19| . 
gives sufficient conditions ensuring almost global stability of an equilibrium point 
for a C^—vector held in R” {i.e., all trajectories, except for a set of initial states 
with zero Lebesgue measure, converge to the equilibrium point). Recently, this tool 
was used by R. Potrie and P. Monzon m to construct a vector held X in R^ where 
the origin is almost globally stable but is not a local attractor for the differential 
system generated for X. We point out that X is an almost Hurwitz held (Hurwitz 
vector held except in a zero Lebesgue measure set). 

This article is focused on two tasks: Firstly, we will construct (a formal de¬ 
scription will be given later) polynomial maps F = XI + H: R^ with JH 

nilpotent, such that the WMYC and the Jacobian Conjecture are true, giving 
the inverse of F explicitly. The results obtained are strongly related with the works 
[4] of L.A. Campbell and |6] of M. Chamberland and A. van den Essen. Secondly, 
we construct two families of three dimensional vector helds having the Rantzer’s 
density functions stated above. The vector helds of the hrst family are a general¬ 
ization of the Potrie-Monzon’s example [18] in the sense that are almost Hurwitz 
and the vector held restricted to the invariant plane z = 0 is a centre. Moreover, 
perturbing these vector helds by XI, we obtain a new family a Hurwitz vector helds 
with the origin as global attractor which are not included in the examples of [121 
Theorem 2.7] and [51 Theorem 2.5]. The vector helds of the second family are not 
almost Hurwitz and the the existence of a density function seems to be the only 
way to demonstrate the almost global stability of the origin. 

The paper is organized as follows. The Section 2 is devoted to generalize the 
procedure of L.A. Campbell in order to obtain new examples to WMYC and the 
Jacobian Conjecture on R”. The Section 3, using density functions, shows two 
family of vector helds such that the origin is almost globally stable. Moreover, we 
construct new examples to MYC. 

2. Examples to Jacobian conjecture and WMYC 

In dimension three, all existing examples and counterexamples to MYC are maps 
of the form XI + H with A < 0 and JH nilpotent (see [8],[7], [5]). This kind of 
maps, in any dimension, also are important in the study of the Jacobian conjecture 
since H. Bass, E. Connell and D. Wright showed in [I] that it suffices to solve 
this conjecture for such maps. Indeed, the Jacobian conjecture follows from the 
injectivity of these maps for all dimensions. 

It is worth to emphasize that the result above triggered new questions and prob¬ 
lems as the following one. Let k be a held of characteristic zero. 

(Homogeneous) dependence problem. Let H = {Hi, ..., i7„) G k \ xi , ... , x^,]" 
(homogeneous of degree d > I) such that JH is nilpotent and H{0) = 0. Does it 
follow that Hi, ..., Hn are linearly dependent over k? 

The homogeneous problem is true in dimension three ([2]) and false in dimen¬ 
sions bigger than hve m)- A counterexample for the inhomogeneous problem in 
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dimension three was given by E. Hubbers in [9]. The map 

H = {y-x^,z + 2x{y - x"^),-{y - x^f ) 

verifies that JJI is nilpotent, rank( Ji?) = 2 and Hi,H 2 , are linearly independent 
over K. Moreover, L.A. Campbell in [4] generalizes this counterexample obtaining 

H = {(t){y - x^), z + 2x(j){y - a;^), -{(j){y - x‘^)f) 

with the same properties. Here (f> is a. C^— function of a single variable. Notice 
that the inverse ot F = I + H can be computed explicitly as follows 

F~'^ = {x-4>{y-x^-z),y-z-2x4>{y-x^-z) + {4>{y-x^-z)f,z+{(t>{y-x^-z)f ). 

Therefore F = I + H with (/) G ]R[t] is an example in dimension three to 

Jacobian Conjecture on R”. Every polynomial map F : M" —)• R" such that 
det JF = 1 is a bijective map with a polynomial inverse. 

Let us emphasize that the previous examples have the special form 
H = {u{x, y, z),v(x, y, z),h{u, v)), 

which is completely studied in [^, which deals with polynomial maps satisfying 
H{0) = 0, h has no linear part and the components of iJ are linearly independent 
over K. Given A = VxUz — UxVz and B = VyUz — UyVz, they show that if JFI is 
nilpotent and deg^ uA ^ deg^ vB, then there exists T G GLz{k) such that TF[T~^ 
takes the form 

( 1 ) {g{t), viz - {bi + 2viax)g{t),a{g{t))'^), 

with t = y + bix + viax"^, via ^ 0 and g{t) G k[t],g{0) = 0 and deg^ g{t) > 1. 

The expressions A and B are very useful for determine if a map F[ is nilpotent 
[6l Proposition 3.1]. Also they are useful for decide if the rows of H are linearly 
independent over k. In fact, we obtain the following result with n = R. 

Proposition 2.1. Let FI = {u(x,y, z),v(x,y, z),h{u{x,y, z),v{x,y, z))) be a poly¬ 
nomial map sueh that the components of FI are linearly dependent over R and h 
has no linear part in u and v. Then A = B = D. 

Proof. Let a, /I G R* such that 

bu * {nxt Jly? by ' {Vx^ Vyj Vz') = O • {Uxi T /? * (Ua,, Uy, Vz^ 
which is equivalent to 


(2) {by o) * {Ux^Uy^Uz) F {by /3) * {Vx^Vy^Vz') - 0. 

Thus, we have the following systems of equations 

{by - a)ux + {by - I3)vx = 

{by - a)Uy + {by - /3)Vy = 

{by - a)uz + {by - I3)vz = 

By using (l3|), we see that 


( 3 ) 


0 

0 

0 . 
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UzA - VzB 


ulvx - UzVziUx + Vy) + vluy 


+ [hy - P){hu - o){ux + Vy) + {hu - a)'^Uy} 

- P)'^Vx 


-{hy - P){hu - a){ ^lllJj) Vx - Vy) - '^\ y} 


(h^-a) 


= 0 


over the set of points where /i„ ^ a. Otherwise, ii hy = a then in ([2]) we have 
that (vx, Vy, Vz) = 0 or hy = p. The first case immediately implies that A = B = 0, 
the second case is not possible due to is a contradiction with no linearity of h with 
respect to u and v. 

In similar way is proved that UzA + VzB = 0 and the result follows. 

□ 


Theorem 2.2. Consider a polynomial map of the form H = (u,v, h(u,v)) such 
that i^(0) =0, h has no linear part and the components of H are linearly indepen¬ 
dent over k. Then if JH is nilpotent and deg^ uA deg,^ vB, for all \ Q the 
polynomial map XI H is injective and has inverse polynomial. 

Proof. According to [6l Corollary 4.1], we can suppose that the components Hi, H2, H3 
of H are as in O- Namely, 


( 4 ) 


If 


Hi{x,y,z) 

H2{x,y,z) 

H'i{x,y,z) 


g{y + bix + viax ^), 

viz — {bi + 2viax)g{y + 61 a; + viax ^), 

a{g{y + bix + viax'^))'^. 


ui = Xxg{ybixviax^) 

U2 = Xy I-viz — {bi2viax)g{y I-biXviax^) 
U 3 = Xza{g{ybixViax^))^, 

it is easy to obtain 


( 5 ) 7U2 + 'ybiUi + 'y'^viaul — 'y'^viu^ = y + bix + viax^, 

where 7 = y- Now, put $ = g{y-\-bix-\-viax^) = g{'yu 2 -\-jbiUi-\-^'^viaul — ^'^viU 3 ) 
and observe that XI H = {Xx + $, Ay + viz — (&i + 2viax)^, Xz + By using 
®, we obtain the inverse of this map which is 

( 7 / + P) = ( 7 a; + Pi, 7 y + P 2 , jz + P 3 ) 


where 






INJECTIVITY AND ALMOST GLOBAL STABILITY 


5 


Piix,y,z) 

P 2 {x,y,z) 

( 6 ) 

P 3 {x,y,z) 


~19i.iy + l^ix + 

—^‘^viz + 7(61 + 2 jviax)g{'jy + 761X + j'^viax'^ — j'^viz) 
—'y'^via{g{'jy + 761a; + j'^viax'^ — j'^viz))'^ 

—aj{g{'yy + jbix + j'^viax'^ — j'^viz))'^. 


□ 

Remark 2.3. For A < 0 (resp. X = 1), the polynomial map F = XI + H of the 
above Theorem is an example of the WMYC (resp. the Jacobian conjecture) onM^. 

Remark 2.4. More examples for the both conjectures in dimension n > 4 can be 
constructed consider the following maps of [9l Proposition 7.1.9]; 


i7i(a:i,... ,a;„) = g{x 2 - a{xi)), 

( — 1 )^ 

Hi{xi,...,Xn) = Xi+i + - -- a(a;i))*“L */2 < I < n - 1 , 

{* - 1 )! 

Hn{xi,... ,Xn) = 7 -— - a(a;i))"“^ 

(n- 1 )! 

where a(a;i) G K.[a;i] with deg a = n — 1 and g{t) G k[t\,g{0) = 0 and deg^ g(t) > 1 . 
Following the lines of the proof of the Theorem we must consider for each 

fixed n > A,Ui = Xxi + Hi, i = 1 ,... ,n and 

1 11 ^1 

^ = g{x 2 - a{xi)) =g(^jU 2 - a(ju) - ^ U 3 + (- 1 )” Mj) 

i=4 

for obtain the inverse (polynomial) of XI + H. In fact, the inverse is 7 / + P where 


Pi{xi,. ..,Xn) 


-l^{xi, ■.■,Xn) 


Pi{xi,...,Xn) = --i{xi+i- ^)(7(a;i- 4 >(a;i,...,a;„)))($(a;i,...,a;„))* ^), 

(*- 1 )! 

{' — 1 

P„(a;i,... ,a;„) = -7-4-— a^”“^)(7(a:i - 4 >(a;i,...,a;n))($(a;i,..., 

(n- 1)! 

with 2 < i < n — 1. Therefore, for A < 0 (resp. X = 1), the polynomial map 
F = XI + H is an example of the WMYC (resp. the Jacobian conjecture) on K." 
with n > 4. 


Example 2.5. We consider n = 4, the inverse of the map XI + {Hi, H 2 , H 3 , H 4 ) 
is 7 / + (Pi, P 2 , P 3 , P 4 ) where 
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Pi{xi,X2,X3,Xa) 

P2 {xi,X2,X^,Xa) 

P3,{xi,X2jX3„Xi) 


-7$, 

-7^X3 + 7^x4 - 7(ai + 7202X1 + ^‘^‘ ia 3 x \)^ + 7^(02 + 803X1)$^, 
-7^x4 + 7(02 + 73a3Xi)<i)^ - 7203$^, 


P 4 (xi,X 2 ,X 3 ,X 4 ) = - 703 $^, 

where $ = $(xi,X 2 ,X 3 ,X 4 ) = ( 7 x 2 — 7 x 1 — 7 ^X 3 — 7 ^X 3 — 7 ^X 3 + 7 ^X 4 ). 

Remark 2.6. A example to WMYC and Jacobian Conjecture, in dimension 4, 
which does not belong to the above family of maps is 

F{x, y, z, w) = (Ax + y,Xy + x^ — w, Xz + y'^, Xw + 2y — z). 

Is easy to see that JH is nilpotent and the rows of JH are linearly independent 
over R and has inverse 

P“^(x, y, z, w) = ( 7 (x - (j), (j), 7(0 - 7 (u> - 27 (x - + 7(2 - (jf^))), 

where (j) = "fy — 7^x^ + 7^2 + 7^u> and 7 = 1 /A. Therefore, for A < 0 (resp. X = 1), 
the polynomial map F is an example of the WMYC (resp. the Jacobian conjecture) 
on 

Remark 2.7. In [SI Corollary 4 . 2 ] it is stated that a complete study of maps 
H = (u,v,h{u,v)) with JH is nilpotent will be achieved when considering the case 
degz(uA) = degz{vB). In this context, the following result gives a partial progress to 
this study since we introduce a large family of maps satisfying degz{uA) = degz{vB), 
which through a linear of change of coordinates have the form described in m- 

Proposition 2 . 8 . The map jI+P with P = {Pi, P 2 , P 3 ) as in has the following 
properties: 

(PO) P{x,y,z) has the form {u{x,y, z),v{x,y, z),h{u{x,y, z),v{x,y, z))). 

(PI) The Jacobian matrix JP is nilpotent and their rows are linearly independent 
over R. 

(P 2 ) degz{uA) = degz{vB). 

(P 3 ) Under the linear change of coordinates {u, v, w) = {x,y — ^viz, z) the map 
7/ + P is transformed into 'yl + P where P has the form (u, v, h{u, v)) and 
degz{uA) ^ degz{vB). 

Proof. (PO) We have P = {u, v, h{u, v)) with h{u, v) = . 

(PI) By using an algebraic manipulator we see that (JP)^ = 0 . We have 
P( 0 ) = 0 and P = (Pi, P2,/i(Pi, P2)) with h{Pi,P2) = —^Pf. Moreover, if we 
consider uj{x) = 61 + 27x1 ax, we have 


B 


-i^viuj{x){g'{f)f - 2-i^vlag{t){g'{t)f --f^vig'{t)+ 
-'y^viuj{x){g' {f)f + 2'y^vlag{f){g'{t)f 
-l^vig'{t) 7 ^ 0 , 
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and 

A = -j^vig'(t)(uj(x) - 2jviag{t)) ^ 0, 
and the result follows from Proposition 12.II 
(P2) We have 

u{x, y,z)-A = -2j^vla{g{t))‘^g'{t) + -f^viuj{x)g{t)g'(t) 

and 

u(x, y,z)-B = -f^vlzg'it) - "f^ViU}{x)g{t)g'(t) + -t^vla{g{t)fg'{t). 

Finally, deg-,(iiA) = max{3fc — 1, 2fc — 1} = 3fc — 1 and deg-,(z;i?) = max{fc, 2k — 
1,3/c- 1} = 3/c- 1. 

(P3) By considering the linear change of coordinate proposed, it is easy to see that 
the coordinates Pi, P2, P3 of P are 


Pi{x,y,z) = g{jy+ -jbix+ -f'^viax‘^) 

P2{x,y,z) = + 7(61 + 27x1 ax)g(7?/+ 761X + 7^X1 ax^) 

P3{x,y,z) = -a7(g(7?/ + 76 ix + 7^uiax^))^. 

This coordinates Pi, P 2 , P 3 are the same that m with 7 = 1, thus degz{uA) ^ 
degzivB). 

□ 

Now we consider, in dimension three, polynomial maps of the form F = XI + H 
with H{0) = 0 and JH nilpotent such that the components of H are linearly 
dependent over K. As it was shown in 0 Proposition 2.1], for such vector field 
F = XI -\- H, there exists T e GlsiM.) such that T~^FT = A/ + (P, Q, 0), where 

P{x,y,z) = -b{z) f{a{z)x + b{z)y) + c{z) and 

(7) Q{x,y,z) = a{z) f {a{z) X F b{z) y) + d{z) 

with a,b,c,d G ]R[z] and / G ]R[z][t]. Also it was pointed out in [SJ Proposition 
2.3] that if A ^ 0, this polynomial map is injective and thus satisfies the conclusion 
of WMYC. In fact, if we consider 7 = y, we can find explicitly the inverse of 
XI + (P, Q, 0), with P and Q as in d?]). Since / G ]R[z][t], in what follows we replace 
f{t) by f{t,z). 

Proposition 2.9. For X ^ {), the inverse of XI + (P, (3,0), with P and Q as in 
is 7 / + {R, S, 0), where 7 = y and 

R = ■jb^'jz) f(^"f{a{"fz) X + b{"fz) y — a{'yz)c{'yz) — b{'yz)d{jz)},"fz'^ — "fc{'jz) 

S = —"fa{'jz) f (^j{a{"/z) X + b{-yz) y — afjz)c{'yz) — b{'jz)dfjz)},"/z'^ — "fdfjz). 
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Proof. Putting 

u = Xx — b{z) f{a{z) X + b{z) y, z) + c{z) 

V = Xy + a{z) f{a{z) X + b(z) y, z) + d{z) 
w = Xz, 

we obtain 

(8) a(z)u + b{z)v = X{a{z)x + b{z)y) + a(z)c(z) + b[z)d{z ), 

which implies 

a{z)x + h{z)y = ^{a{z)u + h{z)v — a(z)c(z) — b{z)d{z)y 

Therefore, by consider that z = yic, we can deduce 

u = Xx — b{'yw) f{m{u,v,w),^w) + c{^w) 

V = Xy + 0(^10) f{m{u,v,w),'jw) + d^-jw), 
with m{u,v,w) = j[a{'jw)u + b{jw)v — a{'yw)c{'jw) — b('jw)d{'jw)}. Finally, 

X = 'y{u + b{'yw) f{m{u,v,w),jw) — c{'yw)) 
y = 'y{v — a{'yw) f{m{u,v,w),jw) — d{'yw)). 


□ 

Remark 2.10. Let F = XI+H ^ be a polynomial map verifying H[0) = 0, 
JH nilpotent and their components are linearly dependent over R. If X < 0 (resp. 
X = 1), then the map F is an example to WMYC (resp. the Jacobian conjecture). 

3. Examples of almost global stability 

In Theorem 3.2,Theorem 3.5] we prove that the vector fields in of the form 
F = XI + H, with A < 0, R as in 0 and g{t) = t + A 2 are counterexamples 
to MYC since they have unbounded orbits. Moreover, by following the respective 
proofs, we can deduce the existence of an open set of initial states whose trajectories 
do not converge to the origin. Thus the origin for these vector fields is not almost 
global attractor in the following sense: 

Definition 3.1. Consider the differential equation 
(9) X = F{x) 

where F : R" —>■ R" is a C^ — map and F{0) = 0. We say the origin is an almost 
global attractor if all trajectories, except for a set of initial states with zero Lebesgue 
measure, converge to the origin. 

From Definition 13.11 it arises the following question: 

Question 1. Do there exist counterexamples to the MYC with the origin almost 
global attractor ? 

In [19] . A. Rantzer introduces a new tool, namely, the density functions in order 
to obtain sufficient conditions for almost global stability of an equilibrium point for 
a vector field in R”. 
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Definition 3.2. A density function of ^ is a map p: M" \ {0} —>• [0,+oo), 
integrable outside a ball centered at the origin that satisfies 

[V • pF]{x) > 0 

almost everywhere with respect to K”, where 

V • [pF] = Vp-F + p[V -i^], 

and S/p, V ■ F denote respectively the gradient of p and the divergence of F. 

The main result of A. Rantzer [TH] is the following. 

Theorem 3.3. Given the differential system 

X = F{x), 

where F G C^, F{0) = 0, suppose there exists a density function p : M" \ {0} —>• 
[0, +oo) such that p{x)F(x)/ || x\\ is integrable on {x G K" : || a:|| > 1}. Then, almost 
all trajectories converge to the origin, i.e., the origin is almost globally stable. 

In this context, the next question arises naturally: 

Question 2. Do there exist counterexamples to the MYC that support density 
functions ? 

In this paper the best answer we give is the following family of almost Hurwitz 
vector fields; i.e. of vector fields for which the Hurwitz condition hold over R" — A 
with A a zero Lebesgue measure set. 

Theorem 3.4. Consider the real numbers c < a < 0, b G M., k > 1, and the 

k 

polynomial R{z) = ^ 02 ^- 2 ^* with a 2 i > 0 for i = 1,... ,k. Then 

i=l 

(10) F{x, y, z) = {y, -x, 0) + {{ax + by)R{z), {-bx + cy)R{z), -zR{z)) 

is an almost Hurwitz vector field. Moreover, this vector field has associated the 
density function p{x,y, z) = {x'^+y'^+R{z))~°‘ with a > max{2, 

Proof. The Jacobian matrix of A in a point {x,y, z) is 

/ aR{z) 1 + bR{z) * 

JF{x,y, z) = I —{l + bR{z)) cR{z) * 

\ 0 0 -{R{z) +zR'{z)) 

Then A 3 = —{R{z) + zR'{z)) is an eigenvalue. The others two eigenvalues Ai, A 2 
verify 

Ai + A 2 = (a + c) R{z) and A 1 A 2 = acR{z)^ + (1 + i>R{z)Y ■ 

Therefore, for z 0 (resp. z = 0), we have A 3 < 0 (resp. A 3 = 0) and Ai and A 2 
have negative real part (resp. null real part). In addition, F verifies the Hurwitz 
condition except in the invariant plane z = 0 . 

In what follows we prove that the map 

p{x, y, z) = {x"^ + + R{z ))~°‘, 

under the conditions of the theorem is a density function for the vector field F. 
The condition a > 2 ensures the integrability of p{x, y, z) outside the ball centered 
at the origin of radius one. 
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It remains to prove that V • (pF){x, y, z) is positive almost everywhere in We 
have 

—a 


^Pix,y,z) = 


{2x,2y,R'{z)), 


and 


Then 


+2/2 + i?(z))“+i 
[V • F]{x, 2 /, z) = (a + c - l)R{z) - zR'{z). 


[V ■ pF]{x,y,z) = {Vp ■ F){x,y,z) + p{x,y,z) [V ■ F]{x,y,z) 

— 0 :R(z) ro ('„^2 I 2 \ 


(a:2 + 2/^ + i?(z))“+i 
1 

^ (x2 + 2/2 + i?(z))“ 

1 


[2(aa;2 + cy^) — z R'{z)] 

[(a + c — 1) R{z) — z R'{z)] 


- (a ;2 + y 2 + i?(^)yi+T + cy‘^)R{z) + azR{z)R'{z) 

+ + 2 /^ + R{z)) [(a + c - 1) R{z) - z R'{z)]] 

^ (a ;2 + 2/2 + i?(z))^+i " l)zi?'(z))i?(z) 

+ ((a + c — 1 — 2aa)R{z) — zR'{z))x^ 

+ ((a + c — 1 — 2ac)R{z) — zR'{z))y'^] . 

k k 

Since R{z) = ^ 2 iz‘^^ and zR'{z) = ^ 2 zo 2 iZ^* with a 2 i > 0 for z = 1,..., fc, we 

i—1 i—1 

obtain [V ■ F]{x, z/, z) > 0 for z 7 ^ 0 and [V ■ F]{x, y,z) = 0 for z = 0, if 
a + c — 1 + 2 (q; — 1 )z> 0 , a + c — 1 — 2 aa — 2 z> 0 , o + c —1 — 2 ac — 2 z > 0 , 
for z = 1,..., fc. Then, the proof is finished due to these inequalities are consequence 


of our hypothesis a > max{ 2 , 


g+c—1—2fc 3—g—c 
2g 


}■ 


□ 


This family is a generalization of the following example of R. Potrie and P. 
Monzon f|18)l 


(11) F{x, y, z) = {y- 2xz'^, -x - 2yz^, -z^), 

which has density function p{x, 2/, z) = {x'^ + + z 2 )“'^. 

Corollary 3.5. The vector field F given by m under the conditions of Theorem 
\3.4\ has the origin as an almost global attractor which is not locally asymptotically 
stable. 

Proof. We have F{x, y, 0) = {y, —x, 0), then the origin is not locally asymptotically 
stable. On the other hand, to prove that the origin is almost global attractor we 
use Rantzer’s result fTheorem l3.3l) . Then it is sufficient to show that the condition 
a > 2 ensures the integrability of p{x,y,z)F{x,y,z)/ || {x,y,z)\\ outside the ball 
centered at the origin of radius one. In fact, if we consider fco = max{l,ai}, 
ki = max{l, b^ + + \b\{a — c)} and ^2 = 2 | 6 | + a — c, we have that 
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||F(a;,y,z)f = + R{zf [{a^ + + c^)y'^ + z'^ + 2b{a - c)xy] 

+2R{z) [(a - c)xy + b{x^ + y"^)] 

< x^ + y^ + R{z f [{b'^ + c^){x^ + 2 /^) + + I b\ {a - c){x^ + y^)] 

+i?(z)(a — c + 2 I 5|)(a;^ + 

< {^2 + y2 + ^2) ^ ^ ^2 ^ i?(z))2(a:" + y^ + ^2^ 

+ fc 2 (:r 2 +y 2 + i?(z))(x 2 + y 2 + z 2 ), 

and that + y^ + i?(z) > feg- Thus, this facts combined with the assumption over 
a imply that 

II f'(x,y,z)f p(x,y,zf _ 1 __fci_ 

||(a;,y,z)f “ (x^ + y 2 + i?(z)) 2 “ (x^ + y 2 + i?(z)) 2“-2 

+ _ ^ _ 

(x 2 + y 2 + i?(z)) 2 “-l 

Therefore p{x,y, z)F{x,y, z)/ || (x,y, z)|| is integrable outside the ball centered at 
the origin of radius one. □ 

Remark 3.6. An alternative -and very simple- proof of Corollary \S.5\ is obtained 
by considering the fact that (F(x, y, z), (x, y, z)) = (ax^ + 6 y^ — z'^)R(z) < 0 for all 
z 7 ^ 0 combined with the invariance of the plane z = 0. 

Motivated by the last remark, we have the following result. 

Proposition 3.7. Let F : R.^ —>• be a C^ —vector field with F(0) = 0 such that 

the plane z = 0 is invariant and [V ■ F](x,y, z) < 0 for all z ^ 0. Suppose that 
there exists a positive C^—function p : R.^ — {0} —>• R satisfying: 

a) [V • pF] (x, y, z) > 0 for all z ^ 0, 

b) limp^o p{p) = oo. 

c) lim||p||^oo p{p) = 0 . 

Then, for any initial state a(0) = (x(0), y(0), z(0)) verifying z(0) ^ 0, the trajectory 
aft) exists for all t G [ 0 , oo[ and tends to zero as t ^ oo. 

Proof. Condition a) implies [V ■ pF](x, y,z) > 0 everywhere, and (Vp ■ F){x, y, z) > 
—p(x, y, z) [V ■ F](x, y, z) > 0 for all (x, y, z) with z 7 ^ 0. 

Let us consider the function V : R^ —?> R defined by C = P~^ outside the origin 
and C(0) = 0. Then V is continuous by b), it is outside the origin and 

1) V(p) >0 for all p and V(p) = 0 if and only if p = 0. 

2) (VV ■ F)(x, y, z) = -p(x, y, z)"^ (Vp ■ F)(x, y, z) < 0, for z 7 ^ 0. 

3) lim||p||_,.oo V(p) = 00 . 

Then, if a(0) = (x(0), y(0), z(0)) is a initial state with z(0) 7 ^ 0, we have -^Viaft)) < 
0 by 2). Then the trajectory a{t) remains over the sublevel C(a(0)), which is 
bounded by 3). Therefore, a{t) exists for all t G [0, oo[ and tend to 0 as t ^ 00 . □ 

Corollary 3.8. If in Proposition \3.T\ we have [V ■ pF](x,y,z) > 0 for all z S R, 
then the origin is globally asymptotically stable. 

If we perturb the vector field F defined by cni) with XI and A < 0, then G = 
A/+F is a Hurwitz vector field, which is another example of MYC, by the following 
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Proposition 3.9. Let G = XI + F with A < 0 and F as in m that verify the 
conditions of Theorem \3.4\ Then G is Hurwitz, admits the same density function 
p{x, y, z) = (a;^+y^+i?(z))“ but with a > 3, and the origin is globally asymptotically 
stable. 

Proof. G is Hurwitz due to the fact that F is almost Hurwitz and JG = XI + JF. 
Since 

\ r rnu \ . {5 — 2a){x‘^ + y"^) + 3R{z) — azR'{z) 

[V-pG]{x,y,z) = [V-pF]{x,y,z) + X^- - %2 ^ y2 + ^ ’ 

under the additional condition a > 3, p is a density function for G. Finally, the 
global asymptotic stability of the origin follows from Corollarv l3.8l □ 

The previous results imply that if we are interested in finding vector fields F 
with the plane z = 0 invariant supporting density functions (but without obvious 
Lyapunov functions), these fields must verify V • -F > 0 over some open set. The 
following family of vector fields satisfy these conditions. 

Proposition 3.10. Consider the vector field F = {P,Q,R) defined by 

P{x, y, z) = —X + H 2 y + oi + 3 02 y + 3 03 a: + 04 + 3 05 z + 

6 06 a: y 2 : + 3 07 y^ z + 3 Os X + 3 og y + oio z^ , 

—A 2 X — y + 61 + 3 &2 y + 3 63 X y^ + 64 y^ + 3 &5 z + 

6 65 X y z + 3 67 y^ z + 3 &8 X z^ + 3 69 y z^ + 610 z^ , 
z (—1 + 3 C 5 x^ + 6 C 6 X y + 3 C 7 y^ + 3 cs X z + 3 cg y z + cio z^), 


Q{x,y, 


R{x,y,z) 

with 


2 3 

08 = 2 — 03 , Oio = 2 (“5 — O7) 


bi 

be 

b9 

C 5 


= —3 02 + 04 + 3 &3 , &2 = 1 — 03 — 


1 


Cio 


64 — 3 — Ol — Cio : 


= o “ 0,7) , = 2 06 + t»5 , &8 = 2 02 — O4 — Og — 63 , 


ClO 


— ~ g + 03 H- — ; ^10 — 3 06 , 

2 Ol 


Cio 


Cio 2 o 4 4 Ol 

03 + , C6 = — - 02 + b3 , C7 = - - y - 03 - ^ 


1 


C8 — — (05 + O7) , Cg — 06 + 65 . 

Then F has the density function 

p(x,y,z) = 


1 


(x^ + y^ + z'^Y 

Moreover the plane z = 0 is invariant and the coefficients of F can be chosen such 
that V ■ F is positive over some open set. 

Proof. Clearly the plane z = 0 is invariant. After straightforward computations we 
obtain 

1 + x^ + y^ + z^ 

(x^ + y^ + z'^Y ' 


[V • pF](x,y,z) 
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and that the integrability condition of p holds. Moreover 

V • F (x, y, z) = —3 + (1 + 4ai)x^ + (13 — 4ai — 4cio)i/^ + + 12b^)yz 

+ (1 + 4cio)2;^ + x((4a4 + I 263 )?/ + 12 a 5 z), 

assume positive values. For example, if y = 0, we have V • F (x, 0, z) = —3 + (1 + 
4ai)x^ + (1 + 4cio)2:^ + 1205 x 2 : which is positive for x, z sufficient large and Oi > 0 
and 05 < 0. □ 


Remark 3.11. The coefficients of the vector field F = {P,Q^R) of Proposition 
\S.1(]\ were obtained from the equation 


[V ■ pF]{x,y,z) 


1 + x^ + 

(x^ + y^ + z^y 


which is equivalent to 


-4 (xP + yQ + zR) + (x^ + [V • F] = (1 + x^ + + z^) (x^ + + z^). 
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